arXiv:l502.06468v4 [math.AP] 14 Apr 2016 


SOME OBSERVATIONS ON THE GREEN FUNCTION FOR THE 
BALL IN THE FRACTIONAL LAPLACE FRAMEWORK 


CLAUDIA BUCUR 


Claudia Bucur 

Dipartimento di Matematica “Federigo Enriques” 
Universita degli Studi di Milano 
Via Cesare Saldini, 50, 1-20133, Milano, Italy 


Abstract. We consider a fractional Laplace equation and we give a self- 
contained elementary exposition of the representation formula for the Green 
function on the ball. In this exposition, only elementary calculus techniques 
will be used, in particular, no probabilistic methods or computer assisted al¬ 
gebraic manipulations are needed. The main result in itself is not new (see for 
instance [2, 9]), however we believe that the exposition is original and easy to 
follow, hence we hope that this paper will be accessible to a wide audience of 
young researchers and graduate students that want to approach the subject, 
and even to professors that would like to present a complete proof in a PhD 
or Master Degree course. 
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1. Introduction 


The Green function for the ball in a fractional Laplace framework naturally arises 
in the study of the representation formulas for the fractional Laplace equations. In 
particular, in analogy to the classical case of the Laplacian, given an equation 
with a known forcing term on the ball and vanishing Dirichlet data outside the 
ball, the representation formula for the solution is precisely the convolution of the 
Green function with the forcing term. As in the classical case, the Green function 
is introduced in terms of the Poisson kernel. For this, we will provide both the 
representation formulas for the problems 


( 1 . 1 ) 


{—Ayu = 0 in Br, 

u = g in 


and 

( 1 . 2 ) 


(-A)®u = g in B^, 

u = 0 in 


in terms of the fractional Poisson kernel and respectively the Green function. More¬ 
over, we will prove an explicit formula for the Green function on the ball. 

Here follow some notations and a few preliminary notions on the fractional 
Laplace operator and on the four kernels that play particular roles in our study: 
the s-mean kernel, the fundamental solution, the Poisson kernel and the Green 
function. 

We briefly introduce the Schwartz space (refer to [11] for details) as the functional 
space 


5(K") := i/eC'®(M”) s.t. VQ:,/3eN”, sup \xi^D°‘f{x)\ < -too 

( £ceR" 

In other words, the Schwartz space consists of smooth functions whose derivatives 
(including the function itself) decay at infinity faster than any power of x; we say, 
for short, that Schwartz functions are rapidly decreasing. Endowed with the family 
of seminorms 


(1.3) 


[/] 


5(R") = SUP(1 + |x|)^ 


a\^N 


the Schwartz space is a locally convex topological space. We denote by 5'(K") the 
space of tempered distributions, the topological dual of 5(IR"). 

We set the following notations for the Fourier and the inverse Fourier transform 
(see for instance, [8] for details) 

/(a:)e-2-“«dx 

jRri 

respectively 

f{x) = B-^f{x)= \ /(e)e2-“«de 

Jr’i 

Here the original space variable is denoted by x e M" and the frequency variable by 
^ e M”. We recall that the Fourier and the inverse transform are well defined for 
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/ e whereas/(a;) = ^ f){x) =T ^ (j^/)(a:) almost everywhere if both 

/ and / e L^(IR”), and pointwise if / is also continuous. Also for all /, g e 

f m9i0di=\ fiOdiOdC 

Jr" Jr" 

The pointwise product is taken into the convolution product and vice versa, namely 
for all f, ge 

Hf*9)=Hf)H9)- 

On the Schwartz space, the Fourier transform gives a continuous bijection between 
5(]R") and 5(K"). We say that / is the Fourier transform of / in a distributional 
sense, for / that satishes 14 ^^ dx < co for some p e N if, for any p e 5(K”) we 
have that 

(1.4) r f{x)ip{x)dx = r f{x)0{x)dx. 

Jr" Jr" 

We remark that the integral notation is used in a formal manner whenever the 
arguments are not integrable. 

Let s e (0,1) be fixed. We introduce the fractional Laplacian for u belonging to 
the Schwartz space. 


Definition 1.1. The fractional Laplacian of u e is defined as 

u{x) - u{y) 


(1.5) 


{-Ayu{x) := C{n,s)P.V. f 
= C{n^ 5 ) lim 




J 

JR- 


X - 

u{x) - u{y) 

X - p|"+2s 


dy, 




where C{n, s) is a constant depending only on n and s. 


Here, P. V. is a commonly used abbreviation for ” in the principal value sense” 
(as defined by the latter equation). The dimensional constant C{n,s) is given in 
[4], formula 3.2, as : 

( 1 . 6 ) . 

where 771 is, up to rotations, the first coordinate of p e K”. 


Remark 1.2. The definition (1.5) is well posed for smooth functions belonging to 
a weighted space, that we define as follows. For s e (0,1) 

Li(M"):={ueLL(K”) s.t. £ 
endowed naturally with the norm 

IIII - r i^(3^)i 4 

Let e be positive, sufficiently small. Then indeed, for u e Lg(IR") and (or 

(ji,2s+s-i £qj. g ^ 2 ^^ 2 ) in a neighborhood of a; e R", the fractional Laplacian is 
well defined in x as in (1.5). See for the proof Proposition 2.1.4 in [12], where 
an approximation with Schwartz functions is performed. We often use this type 
of regularity to obtain pointwise solutions to the posed equations. We will always 
write (7^®+® to denote both (7°’2®+® for s < 1/2 and for s > 1 / 2 . 
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In definition (1.5), the singular integral can be substituted with a weighted sec¬ 
ond order differential quotient by performing the changes of variables y = x + y 
and y = X — y and summing up. In this way, the singularity at the origin can be 
removed, as we observe in the following lemma (see for the proof Section 3.1 in [3]). 


Lemma 1.3. Let (—A)® be the fractional Laplace operator defined by (1.5). Then 
for any smooth u 


(1.7) 


i-AYnix) 


C{n,s) r 2u{x) 

2 Jb’* 


u{x + y) — u{x 

|y|n+2s 


y) 


dy. 


For u e iS(K") the fractional Laplace operator can also be expressed by means 
of the Fourier transform, according to the following lemma (see Proposition 3.3 in 
[4] for the proof). 


Lemma 1.4. Let (—A)® be the fractional Laplace operator defined by (1.5). Then 
for any u G 5(]R"’) 


( 1 . 8 ) (-A)®u(a:) = J-i((27r|e|)^*S(0). 


We refer usually to pointwise solutions, nevertheless distributional solution will 
also be employed. Following the approach in [12] (see Definition 2.1.3), we introduce 
a suitable functional space where distributional solutions can be defined. Let 

5«(K”) := {/gC'”(M”) s.t. VogN”, sup {l + \xY+'^’^)\D‘^f{x)\ < -hcx)}. 

The linear space 5s (K") endowed with the family of seminorms 

[/]5,(B") := sup (1-h |a;r+^®)|D“/(a;)| 

xeB" 

is a locally convex topological space. We denote with 5((K”) the topological dual 
of 5s(K"). 

We notice that \iip e 5(K") then (—A)®(^ g 5s(ffi"'), which makes this framework 
appropriate for the distributional formulation. In order to prove this, we observe 
that for any x G the bound 

(1.9) |(-A)®(/j(a;)| ^ 


follows from the upcoming computation and the fact that ip G 5(]R") 

l(-A)V(x)| 




I 


\2p{x) - p(x -y)- p{x + y)\ 


\y\ 


n+2s 


dy + 2 


\ 

JR-' 


\ip{x) -ip{x + y)\ 




\y\ 


n+2s 


dy 


^ sup (1 -b \z\r+^\D^p{z)\ + sup (1 + \z\rMz)\ + ||<b>||LqB") 

\ zeR-^ 

Moreover, we observe that, up to constants, 

d„,(-A)®(/j(x) = d,^,J-^(^|^p®(^(0)(a:) = .5"^ (^i^i|Cp®^(C)) (x) 

= = i-AYd^M^). 

Hence, by iterating the presented argument, one proves that (— A)®(/j g 5s(ffi"'), 
which assures our claim. Then we have the following definition. 
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Definition 1.5. Let / e we say that u e 5'(M") is a distributional solution 

of 

(-A)"u = / in M” 
if 

(1.10) < u, (—A)^(p >s= I f{x)(p{x)dx for any e 

Jb" 

where < >s denotes the duality pairing of 5'(K") and 5s(K”) and the latter 
(formal) integral notation designates the pairing 5(K") and 5'(K"). 

We also use the integral notation for the pairing < >s in a purely formal 
manner whenever the arguments are not integrable. We notice that the inclusion 
Lg(K"') c 5'(K"') holds, in particular for any u e Lg(K") and e 5s(K") we have 

( 1 . 11 ) 


< U,tp >s 


< 


r |u(a;)| \'ilj{x)\dx ^ f 
Jr" Jr' 


i{a 


{l+\x\^+^^)\ijix)\dx 


l + |a;|"+2® 

^ [V']5,(R")ll^llil(R")- 

We introduce now the four functions Ar, Pr and G, namely the s-mean kernel, 
the fundamental solution, the Poisson kernel and the Green function. The reader 
can see Section 2.2 in [7] for the theory in the classical case. 

Definition 1.6. Let r > 0 be fixed. The function Ar is defined by 


( 1 . 12 ) 


^r{y) ■■ = 


z{n, s) 


(|i/|2 - r2)®|?/|" 


0 


y e E.'^\Br, 
y G Br^ 


where c(n, s) is a constant depending only on n and s. 
Definition 1.7. For any x e IR"\{0} the function if* is defined by 


(1.13) 


$(x) := 


z(n,s)|; 
1 ' 


I— n+2s 


if n 7 ^ 2s, 


o(l,-)log|x| ifn = 2s, 

where a(n, s) is a constant depending only on n and s. 

Definition 1.8. Let r > 0 be fixed. For any x e Br and any y e 
Poisson kernel Pr is defined by 


^\Br, the 


(1.14) 


Priy,x) := c{n,s) 


— IxP 


1 


\y\2 — j |a; _ y|Tl 


The Poisson kernel is used to build a function which is known outside the ball 
and s-harmonic in Br- Indeed, it is used to give the representation formula for 
problem (1.1), as stated in the following theorem. 

Theorem 2.10. Let r > 0, g e Lg(IR”) n C'(M") and let 

, , Pr{y,x)g{y)dy ifxeBr, 

Ug{x) := < Jr"\_b,, 

[g(x) i/xeIR"\i?r. 


(1.15) 
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Then Ug is the unique pointwise continuous solution of the problem (1.1) 


(—A)®m = 0 in Br, 

u = g in W^\Br. 


Definition 1.9. Let r > 0 be fixed. For any x,z e B^ and x ^ z, the function G 
is defined by 

(1.16) G{x,z) ■.= ^{x - z) - [ ^{z-y)Pr{y,x)dy. 


From this definition, a formula that is more suitable for applications can be 
deduced. Indeed, one of the main goals is to prove this simpler, explicit formula 
for the Green function on the ball, by means of elementary calculus techniques. At 
this purpose, Theorem 3.1 establishes a symmetrical expression for G. 

Theorem 3.1. Letr > 0 be fixed and let G be the function defined in (1.16). Then 
if n 2s 


(1.17) 


G{x, z) = K{n, s)\z — X 


2s—n 


I 


roix,z) fs-l 


(t + 1) : 


dt, 


where 

(1.18) 




and K{n, s) is a constant depending only on n and s. 
For n = 2s, the following holds 


(1.19) 


GO. 0-^(1 ,log - 



This result is not new (see [2, 9]), however, the proof we provide uses only 
calculus techniques, therefore we hope it will be accessible to a wide audience. It 
makes elementary use of special functions like the Euler-Gamma function, the Beta 
and the hypergeometric function, that are introduced in the Appendix (see the 
book [1] for details). Moreover, the point inversion transformations and some basic 
calculus facts, that are also outlined in the Appendix, are used in the course of this 
proof. 

The main property of the Green function is stated in the upcoming Theorem 3.2. 
The function G is used to build the solution of an equation with a given forcing 
term in a ball and vanishing Dirichlet data outside the ball, by convolution with 
the forcing term. While this convolution property in itself may be easily guessed 
from the superposition effect induced by the linear character of the equation, the 
main property that we point out is that the convolution kernel is explicitly given 
by the function G. 

Theorem 3.2. Let r > 0, h e C'^®+®(ilr.) n G{Br) and let 


g , ) Hy)G{x,y)dy ifxeB^, 

u{x) := < Jb,. 

lo 


ifxeBT\Br. 
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Then u is the unique pointwise continuous solution of the problem (1.2) 

J (— A)^m = h in Br, 

= 0 in MT'\Br. 

The proof is classical, and makes use of the properties and representation for¬ 
mulas involving the two functions $ and Pr- 

We are also interested in the values of the normalization constants that appear in 
the definitions of the s-mean kernel (and the Poisson kernel) and of the fundamental 
solution. We will deal separately with the two cases n ^ 2s and n = 2s. We have 
the following definition: 


Definition 1.10. The constant a{n,s) introduced in definition (1.13) is 

r(f-s) 


( 1 . 20 ) 

( 1 . 21 ) 


a(n,s) : = 


22«7rtr(s) 


for n 7 ^ 2s, 




for n = 2s. 


( 1 . 22 ) 


The constant c{n,s) introduced in definition (1.12) is 

r(^) sinTTs 


c(n, s) := 


These constants are used for normalization purposes, and we explicitly clarify 
how their values arise. However, these values are only needed to compute the 
constant /c(n, s) introduced in Theorem 3.1, and have no role for the rest of our 
discussion. The value of K(n, s) is given in terms of the Euler-Gamma function as: 
Theorem 3.3. The constant K{n, s) introduced in identity (1.17) is 

r(-) 

= 02s I ^ 

2^^7r 2 1 ^(s) 

k(i, ^ forn = 2s. 

One interesting thing that we want to point out here is related to the two con¬ 
stants C(ri,s) and c(u, s). The constant C{n,s) is given in [4] in the definition of 
the fractional Laplacian, and we defined it here in (1.6). The constant c(n, s) is 
introduced in [10] in the definition of the s-mean kernel and the Poisson-kernel, 
and is here given in (1.22). These two constants are used for different normaliza¬ 
tion purposes. We observe that they have similar asymptotic properties, however 
they are not equal. In the following proposition we explicitly compute the value of 
C{n,s) defined in (1.6). 

Theorem 3.5. The constant C(n, s) introduced in (1.6) is given by 


(1.23) 


C{n, s) 


2^"sr (f s) 
7rtr(l — s) 


This paper is structured as follows: in the Preliminaries (Section 2) we introduce 
some kernels related to the fractional Laplacian. In Subsection 2.1 we define the 
s-mean value property by means of the s-mean kernel and prove that if a function 
is s-harmonic, then it has the s-mean value property. Subsection 2.2 deals with 
the study of the function $ as the fundamental solution of the fractional Laplacian. 
The fractional Poisson kernel is introduced in Subsection 2.3, and the representation 
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formula for equation (1.1) is obtained. Section 3 focuses on the Green function, 
presenting two main theorems: Theorem 3.1 gives a more basic formula of the 
function G for the ball and is treated in Subsection 3.1; Theorem 3.2, that illustrates 
how the solution to equation (1.2) is built by means of the function G, is dealt with 
in Subsection 3.2. The computation of the normalization constants introduced 
along this notes is done in Subsection 3.3. The Appendix A introduces three special 
functions (Gamma, Beta and hypergeometric), the point inversion transformations 
and some calculus identities that we use throughout this paper. 


2. Preliminaries 

In this section, we deal with the s-mean kernel, the fundamental solution and 
the Poisson kernel. 

Let s e (0,1). 


2.1. The s-mean value property. In this subsection, we define the s-mean value 
property of the function u, namely, an average property defined by convolution of 
u with the s-mean kernel. 


Definition 2.1 (s-mean value property). Let x e M". We say that u belonging to 
Lg(IR"') and continuous in a neighborhood of x has the s-mean value property at x 
if, for any r > 0 arbitrarily small, 

(2.1) u{x) = Ar * u{x). 

We say that u has the s-mean value property in 17 c K" if for any r > 0 arbitrarily 
small, identity (2.1) is satisfied at any point x e fl. 

The kernel defined in (1.12) is used to state the s-mean value property, which 
makes it reasonable to say that A^ plays the role of the s-mean kernel. The main 
result that we state here is that if a function has the s-mean value property, then 
it is s-harmonic (i.e u satisfies the classical relation (—A)®m = 0). 

Theorem 2.2. Let u e L)(]R") be G^®+‘^ in a neighborhood of x e K". If u has the 
s-mean value property at x, then u is s-harmonic at x. 


Proof. The function u has the s-mean value property for any r > 0 arbitrarily 
small, namely 


u{x) = Ar * u{x) = Ar{y)u{x — y) dy. 

Using identity (A. 19) we obtain that 


0 = u{x) - r A,{y)u{x -y)dy = c(n, s)r2* f 
Jb"\b.. Ji 

u{x) — u{x — y) 


thus, since r > 0 

( 2 . 2 ) 


J 

Jr- 


\B^ {\y? -r^Y\y\ 

Hence, in order to obtain (—A)®u(a;) = 0 we prove that 


u{x) — u{x — y) 
\b/(M^ - r2)%|’ 

dy = 0. 


-dy, 


(2.3) 


lim 

r—»-0 


J 

JR- 


-\B, 


u{x)-u{x-y) 

- JPP2S - dy = hm^ 


\y\' 


im r 

^0 Jb' 


u{x) — u{x — y) 

\B, {\y?-r'^Y\yY 


dy. 
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Let R > rVS- We write the integral in (2.2) as 

u{x) - u{x -y) , 
ay 


J 

Jr- 


(2.4) 


\B^ {\y? - r'^Y\y\'^ 

u{x) — u{x — y) 

-\br {\y\^ - r'^YlvY 

= h{r,R) + l 2 {r, R). 


j 

Jr- 


dy + 


I 


Br\b^ 


u{x) - u{x -y) , 

{\y?-rY^\yY ^ 


In Ii{r,R) we see that 


\v\^ 

|y|2_^2 


< 2 and obtain that 


\u{x) - u{x - ^)| 

(|y|2 _ 


< ns Hx) - u{x - y)\ 

'' |y|n+2s 


eLYR^\BR,dy), 


as u e Lg(K"). We can use the dominated convergence theorem, send r ^ 0 and 
conclude that 


(2.5) 


lim/i(r, i?)= I 

’■^0 Jr> 


u{x) — u{x — y) 


■\Br 


\y\ 


n+2s 


dy. 


Now, for r < \y\ < R and u e 0"^^+^ (for s < 1/2) in a neighborhood of x we have 
the bound 


\u{x) - u{x - y)\ < c|j/p®+®, 
while for s > 1/2 and u e ( 7 i. 2 s+e-i ^gg 


{x) — u{x — y) — y ■ \/u{x)\ = y(\^u{x — ty) —'Vu{x)) 

Jo 

^ \y\ Vu{x — ty) — Vu{x) 

Jo 


< c(s,e)|y| 


2s+e 


Notice that (\y\'^"^Y^Y\y\n s-iid are even functions, hence they vanish when 

integrated on the symmetrical domain BR\Br- Therefore, by setting 


( 2 . 6 ) 

we have that 
J(r,i?) = r 

Jl 


J{r,R) :=h{r,R)- f 

Jl 


i{x) - u{x - y) 


dy 


br\b^ |yp*+” 

u{x)—u{x — y) — y ■ Vu(x) u(x)—u(x — y) — y ■ Vu(x) 


Br\B^ 


2 _ J^2 ^ S I rjj I Th 


\y\ 


2s+n 


dy 


and by passing to polar coordinates and afterwards making the change of variables 
p = rt we get 

|J(r,i?)|^c(s,£) r \yr^H{\y?-Y)-^\y\---\y\-^-^^)dy 

JBr\B^ ^ ' 

= c(n,s,e)J^ 

< c{n,s,e)r^ dt 
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since t/{t + 1) > 1. Now for t e (1, V^) we have that 


- l) dt < c(s) (^(t - 1)-* - dt = c{s). 


On the other hand, for t > -\/2 


and we have that 

R 


lim 
1 — 


fj-i 


/ In -s 

s ( 

1 

(1 - - - 

1 < -(1 - 


V tJ 

t V 

V2 


pGO 


\dt^ \ 


(t-iy ) 

Jv2 

1 

pOO 


< c(s) 



JV2 


1 \-«-i 


dt 


Thus by sending r —> 0 we obtain that 

and therefore in (2.6) 


lim J(r, R) = 0 

r—>0 


lim I 2 
1 — 


(r,R) = lim f 


i{x) - u{x - y) 


Br\B^ 


\y\ 


2s+n 


dy. 


Using this and (2.5) and passing to the limit in (2.4), claim (2.3) follows and hence 
the conclusion that (—A)®M(a:) = 0. □ 


2.2. The fundamental solution. We claim that the function $ plays the role of 
the fundamental solution of the fractional Laplacian, namely the fractional Lapla- 
cian of $ is equal in the distributional sense to the Dirac Delta function evaluated 
at zero. The following theorem provides the motivation for this claim. 

Theorem 2.3. In the distributional sense (given by definition (1.10)^ 

(-A)«$ = do. 


The computation of the Fourier transform of the fundamental solution is required 
in order to prove Theorem 2.3. 


Proposition 2.4. a) For n > 2s, let / e L^(IR”) n C(M”) with / e 5s(M”), 

b) for n ^ 2s, let / e L^(IR) n (7(M) n ((— 00 ,0) u (0, + 00 )) with / G 5s(K) such 

that 


(2.7) 


|/(a:)| ^ ci|x|' 

i/wi < g 

l/'MI < c[lx 


2s-l 




for a; G M 
for |x| > 1 
for 0 < |a;| ^ 1 
for |x| > 1. 


Then 


r ^{x)f{x)dx= r (27r|a;|) ^^f{x)dx. 

Jb"- Jb" 
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Proof. We notice that the conditions (2.7) on / assure that the integrals are well 
defined. Indeed, since $ e c the left hand side is finite thanks to 

(1.11). The right hand side is also finite since, for n > 2s, 

f l/(a;)||a;r^®da: ^ c„ sup |/(a;)| r |/(a;)||a:|“^® dx 

Jk" xeBi Jo Jh"\Bi 


xeBi 

^ Cn sup |/(x)| + ||/||l1(E-) 
x^Bi 


and for n ^ 2s we have that 

dx ^ 


r \f{x)\\x\-^^dx^ \ i/c 

Jr Jr\Bi 


I 


a;)||a;| ^ da; + Ci | dx 

Bi 

< II/IIli(®) + ■ 

a) For n > 2s we prove that 

(2.8) a(n, s) f |a;|“"'"'"^®/(a;) dx = f (27r|x|)“^'^/(x) dx. 

We use the Fourier transform of the Gaussian distribution as the starting point of 
the proof. For any d > 0 we have that 

jr(g-^5ft|=) ^ 

In particular for any / e L^(M") and / e 5s(K”) (which is a subspace of L^(IR")), 
by Parseval identity we obtain 

r e-’^'5|^l'/(x)dx = f S-^e-^^f{x)dx. 

JR" JR" 

We multiply both sides by integrate in 5 from 0 to oo. We use the notations 


h = 


J” J d5-®-ie-"^l"l'/(a;) dx^ dS 


and 


h = 


rfj 

Jo \ Jr' 


,5 s ig ■5 /(x) dx dS, 


having Ii = 12- We perform in Ii the change of variable a = djxp and obtain that 

r-OO \ 


h = 


f \xr+^^f{x)(r 

\ JO 


a 2 




We set 
(2.9) 


Cl := 


-f 


a 2 


^ —s — l^—TTa 


da, 


which is a finite quantity since — s — 1 > —1. On the other hand in I 2 we change 
the variable a = |xp/d and obtain that 


I 2 = 


We then set 

( 2 . 10 ) 


r \xr^f{x)(ra^-^ 

JR^ \ Jo 

rCO 

C2 := 

Jo 


da, 
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which is finite since s—1>—1. As/i=/ 2 it yields that 

f |a;r"+^V(a;)da;= r {2Tr\x\)-^^f (x) dx. 
Jr" Jr" 


We take 

( 2 . 11 ) 


a(n, s) = 


Cl 


C2{2tt) 


2s 


and the claim (2.8) follows. This concludes the proof for n > 2s. 

b) For n < 2s (hence n = 1 and s > 1/2), let i? > 0 be as large as we wish (we 
will make i? go to oo in sequel). Then 


L 


,|2s-l 


J{x) dx = ^ ^ (J{x) + f{-x)^ dx 

= 2 r r /(5) cos(27r^a;) d^ dx 

Jo Jr 


= 2 


J fiO ^ J x^^ ^ cos( 27 r^a;) dj^ d^. 


We use the change of variables x = 2 ttx (but still write x as the variable of inte¬ 
gration for simplicity) and let R = 2TrR. Then 

rR rR 

cos(27r^a;) dx = (27r)“^® cos(^x) dx. 

Jo Jo 


We have that 




|xp® V(a:)dx = 


J (/' 


x^® ^ cos(^x) dx ) d^. 


Integrating by parts and changing variables |5|x = t we obtain that 


r-R 

x^® ^ cos(^x) dx 

Jo 


= X 


2s-iSin(Cx) 


„2s-2sin(Ca;) 




-(2s-l) 

JO 

(•m 

(2s-l)|Cr^M t^®"^sintdt. 
Jo 


Therefore 

( 2 . 12 ) 


L 


|x|2®-7(x) dx = 

TJ.2S 


sin(^i?) 


de 


f /(e) , 


m\ 


t^® ^sintdtld^. 


We claim that 
(2.13) 


lim R 
i?—»-00 


2s-l 


f no 

Jr 


sin(^i?) 

e 


d^ = 0. 
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We integrate by parts and obtain that 

J t?, 

Jo 


rCC 

mr- 

JO 




< 


1/(61 

cos(^J?) 


R 

rOD 


60 

cos{^R)\ 


1 ^ 
^ R 


I 

Jo 


cosm\\lMd^ 
0 s 


By (2.7), for ^ large we have that 

hence li„ !M!IffSp! . o. 




R 


For ^ small we have that 
1/(61 


, hence lim^MI^.O. 

^ €-0 ^ i? 

Furthermore, by changing variables t = ^R (and noticing that the constants may 
change value from line to line) we have that 


r 

Jo 


and 


f 


1/(61 , 

6 ' 


If (61 


f' 

Jo 


cos(^-R)| ^ Cl ^1 cos(^i?)| + C 2 ^ '^| cos(^i?)| 


i: 


^ Cii? 


l-2s 

JO 


JR 


^1 cost] dt + C 2 i?^ I t '’|cost|dt^- 


4 f' 

JO 


cos{^R)\^ c[ "^1 cos(^i?)| + C 2 ^ cos(^i?)| d^ 


^ c\R 


6: 

/ d1-2s r t^^-‘^\cost\dt + C 2 R f 

Jn Ji 


t ^1 cost] dt ^ 


Hence 


J 'OO 

/(6 

0 


sin(^i?) 


dC 




and in the same way we obtain 

sin(Cf) 


r = r/(-o 

J -00 s Jo 


i?’ 

sin(^f) 




dC 


< 


R' 


Therefore 


lim i? 2 «- 
R^od 


' f /(6 

M 


sin(^i?) 


dC = 0 


and we have proved the claim (2.13). Now we claim that (and this holds also for 
n = 2s) 

(2.14) 

lim f /( 6 l/r^^( f f^"^sintd6 d^ = -cos(7rs)r(2s-1) r /(^)|^|"^"dC. 

■R-’® Jr ^ Jo '' Jr 

In order to prove this, we estimate the difference 


rCO 

r-RkI 

f^“^sintdt 


^00 

f®“^sintdt 

Jo 

Jo 


Jri«i 


^ |f® ^cost] 




rOO 

( 2 s — 2 ) cost] dt ^ c(f 1 ^ 1 )^®“^. 

Jr|5| 
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We then have that 


shit dt — sint dtj 

»/IR «^0 »-'0 


^ cR 


^ cR^ 


,1s—2 


f \nom-^ 

JR 

Jo Jl 


rfC + C2 r 


Hence we obtain 


^sintdt^ = J^/(^)|^| ^shitd^d^ 

and the claim (2.14) follows from the identity (A.36) in the Appendix. 

By sending R to infinity in (2.12) we finally obtain that 


(2.15) 


r ^/(x)da; = ^-^(2s-l)cos(7rs)r(2s-1) f |^| ^''f{Cld£. 

Jb Jr 

= 2cos(7rs)r(2s) f (2^|e|)-2V(C) d?- 

JR 


Therefore taking a(l,5) = (2cos(7rs)r(2s)) ^ we obtain that 

a(l,s) r \x\^‘'~^J{x)dx = r {2'k\x\)~'^‘' f{x)dx, 

Jr Jr 


hence the result for n < 2s. 

Now, for n = 2s we have that 


I; 


log \x\f{x) dx = 2 


J if 


logaicos(27r^a;) dx d^. 


We change variables x = 2ttx (but still write x as the variable of integration for 
simplicity) and let R = 27ri?. Then we have that 

/ \ dx 

I logxcos(27r^x) dx = I (^logx — log(27r)j cos(^x) — 

271" Jo 


2tt 


c«l?l 

Jo 


sint 


We integrate by parts and obtain that 

r logxcos(^x) dx = — -r^ 

Jo s l?l 

We thus have that 

r iogixi/(x)dx=-logi? r f /(c)icr^ ( f 

JB„ TT Jr 4 TT Jr \ Jq 


dt. 


Sint , \ 

dt j d 4 . 


We claim that 
(2.16) 


lim logi? 

R—>CX) 


[ /(C) 

Jr 


sin(Ci?) 

c 


dC = o. 
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Indeed we have that 

Jo 4 

Moreover integrating by parts we have that 




< 


rl/R 

l/(ai 

Jo 


^ ^ Cii? 


rl/n 

Jo 


^ ^ R 


rOO 

m) 

Jl/R 


sin(^i?) 


dC 


< 


l/(/)l 

COs(Ci?)| 


R 


l/R 




I 


l/R e 


ffl|cos(Ci?)|de 


+ 


/ 


L/M 

l/R i 


cos(^i?)| d^ . 


We have that for ^ large 


1/(01 |cOs(Ci?)| C 2 |cOs(Ci?)| 




R 


< 


R 


e 


hence 


|/( 0 ||cos(ji^)| 

{->■00 ^ R 

On the other hand by using the change of variables t = ^R 


I 

Jl 


l/R / 


I^P|cos(ei?)|de<ci 


I 

r 


|cos(^i?)| 

-a4 + C2 


j: 


C0s(^i?)| 

0 


d^ 


I cost] 


dt + C 2 


J- 

JR 


R^ 


cost] dt ^ Cl logi? + C 2 . 


Moreover 

lf(OI 


/ 


l/R i 


cos(^i?)| d^ ^ d 




l/R C 
^ ' cost] 


cos(ei?)| .. , 

“S + ^2 


r 


C0s(^i?)| 

0 


d^ 


Hence 


lim logi? 
-R—>00 


, r I cost! , , r® i?, ,, 

■l -^ dt + C 2 -^1 cost] dt ^ Cl log i? + C 2 . 

Jl ^ JR i 

rCO 

fiO 

Jo 


sin(^i?) 


dC = 0 


JO / 

and since the same bounds hold for /(O d^, the claim (2.16) follows. 

Also, the proof of the claim (2.14) assures us that 

r/(OIO-'( 

R^JjJo vJo 


It follows that 

hence 

(2.17) 


J 

Jm 

1 

TT 


log|a:|/(a:) dx = 


" r /(OIO”' ( J” ^ 
^00 

= 2^ /(oicr'rfo 

f ierV(Orfo 

Je 


r \og\x\f{x)dx = r (27r|^|) ^fiOd^ 

Jr Jr 


and the result holds for n = 2s. This concludes the proof of the Proposition. □ 
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Remark 2.5. It is now clear that we have chosen a{n, s) in Definition 1.10 in order 
to normalize the Fourier transform of the fundamental solution. Indeed, for n > 2s, 
we perform the change of variable ira = t in (2.9) and by (A.l) we obtain that 


Cl = Tt'* 2 


f 




'e-'dt = 7r"-2r( --S 


Also in (2.10) we change the variable Tra = t and get that 

nOO 

’ dt = '!t~^T{s). 

Jo 


C 2 = TT 


Therefore 


Cl 7r--tr(t-s) r(f-s) 


C2 


r(s) 


hence by (2.11) a{n,s) = 


22*71 2 r(s) 


The value a(l, s) is computed in (2.15). We point out that we can rewrite this value 
using (A.6) and (A.4), as follows 


a(l,s) = 


r(l /2 - s) 


2 cos(7rs)r(2s) 2^®y'^r(s) 

Moreover, we observe that identity (2.17) says that 

By applying this latter Proposition 2.4, we prove Theorem 2.3. 

Proof of Theorem 2.3. For any / e iS(K”) we have that ^|^p®/(^)^ e 5s(K") 

(according to definition (1.8) and to (1.9)). Notice that |^P®/(^) e L^(IR") nC(K"), 
since 

f kPl/(a;)| dx ^ f \x\'^^-'^-^ dx + sup \f{x)\ < c(/), 

Jr" Jr"\Si xeBi 

where we use the seminorm defined in (1.3). Moreover, for n ^ 2s we have that 

l^ri/(ai<ll/llL==(R)ieP^ = ci|eP'^ for^GK, 

\evm ^ [/(e)]5(R)l?P^-" < I for lei > 1. 

Also for 0 7^ ICI < 1 


|(iepv(a)| < 2s\e^-^\m\ + 


< ler-M 2 sii/ii 


\dm 

d^ 


L^( 


= c;iei^ 


and for |^| > 1 


|(iepv(c))| 

< 2s[/]°f^JeP^-^ + < c'lei 
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which proves that / satisfies (2.7). From Proposition 2.4 it follows that 
<$,(-A)V>,= r ^x)T-U{27r\(\f^m){^)dx 


r {27r\^\)-^%2n\^\r^m = f fiO d^ = /(O). 


□ 


Therefore in the distributional sense 

(-A)*<i> = Jo. 

We have the following Lemma. 

Lemma 2.6. Let / e Cc(K”), then / * $ e Lg(IR"). 

Proof. To prove that / *$ e Lg(IR"), we suppose that supp / c Bn and we compute 


J 

JW‘ 


1 + |a;|”+ 2 ® 


dx = 




- y) 

1 + |a;|”+2® 


dx ) dy 


*j Bji \ •JK' 


We set 
(2.18) 


(I. 


and prove it is a finite quantity. We take for simplicity i? = 1 and remark that the 
constants in the next computations may change value from line to line. For n > 2s 
we have that 

Jb. (L ■'") Jb. (L ■'") 


For X small we have that 

|l - yf—' 


f (\ 

JBi \JB2 1 


+ \x 


\n-\-2s 


dx] dy ^ 


Jb.(JbJ 

■"‘Jb. (I 

c„,« r ( 2 + 

do 


yr-^dx] dy 


^ dt] dy 


dt = 


For X large, we use that \x — y\^ |a;| — |y| and 1 + |x|”+^® > thus 


J (( (J, (I 

JBi \Jk"\B 2 + FI / dsi \Jr"\B2 


a;| - |n|)"®-”|x|-"-2® 


dx ) dy 




^ c, 


r dp) 

^00 

,J^ (p-l) 2 ®-"-pp = c„.«. 


dt 
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Hence for n > 2s the quantity Cn^s,B. in (2-18) is finite. Meanwhile, for n < 2s for 
X small the same bound as for n > 2s holds. For x large, we have that 

X — 


I (J 

JBi \JW' 


1 + 


dx ] dy 


2s-li I-1-2S 


{\x\ + \y\r-^\x 


J (I 

JBi \J 1 


dx j dy 
dt = Cs- 


In the case n = 2s from the triangle inequality we have that 

log \x - y\ ^ 

1 + IxP 


and 


J (I 

JBi \JB 

r / r \og\x-y\ 
JBi VJr' 


1 + bl 


dxj dy ^ c J log(t + 1) dt 

rOD 

J log(t + l)t 


= c 


dx] dy ^ 


dt = c. 


Hence Cn,s,R in (2.18) is finite and we have that 


(2.19) 


f l/*^(a;)(x)| 

Jr TTN^ C„.«,«||/Iboo(r). 


It follows that / * $ e Lg(IR"’), as stated. 


□ 


Before continuing with the main result of this section, we introduce the following 
lemma, that will be the main ingredient in the proof of the upcoming Theorem 2.8. 

Lemma 2.7. Let / e (^“(K"), let ip be an arbitrary function such that we have 
p e 5s(K”) and the following hold: 

a) for n > 2s, (/? e L^(IR") n (^(K"), 

b) for n ^ 2s, (fi B L^(IR) n C(M) n ((—QO, 0) u (0, oo))) and 

\p{x)\ < cibp® 


\p{x)\ ^ 


C 2 


\ip\x)\ ^ c\\x 

l‘d'(^)l < A 


2s-l 


for x e K 
for bl > 1 
for 0 < bl ^ 1 
for bl > 1. 


Then 

( 2 . 20 ) 


r f * ^{x)p{x) dx = r (27rb|) f{x)ip{x) dx. 
Jr" Jr" 


Proof. In order to prove identity (2.20) we notice that by the Fubini-Tonelli theorem 
we have that 

/ * ^{x)p{x)dx = ( ^(y)/(ai -y)dy \ p{x) dx 

Hy)i f{x-y)p{x)dx] dy. 

\ Jr^ / 


We denote 


f*p{y)-= f {x - y)p{x) dx = f{x)p{x + y)dx 

Jr^ Jr^ 
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and write 
( 2 . 21 ) 


f * ^{x)ip{x)dx = ^{y)f*0{y)dy. 

Jr" Jr" 


The operation * is well defined for / e C“(K") and ip e 5s(ffi"), furthermore it is 
easy to see that 

J^(f*ip)(x) = f{x)p{x). 

We notice at first that since p and ip are continuous, if (ip) = yi on M". We define 

( 2 . 22 ) tjj{x) := T{fiip){x) = f{x)p{x) 
and we write (2.21) as 

(2.23) r f * ^{x)ip{x)dx = r ^{y)^{y)dy. 

J]R»^ 

To apply Proposition 2.4, we have to check that iIj verifies (2.7). Since 

\^{x)\dx= \ |/(x)||(^(x)| dx ^ ||/||i^(K-)||(/j||Li(R-), 

Jr" Jr" 

we have that ip e L^(IR"). Also, ip e (^(M”) as a product of continuous functions. 
We claim that f*ip e iSs(K"). Indeed, suppose supp/ c Br for i? > 0. We remark 
that in the next computations the constants may change from line to line. Then 
for |x| ^ 2i? we have that 

(1 + |x|”+^^)|/5^(x)| ^ Cn,s,R f |/(y - x)ip{y)\dy 

< Cn,s,R\\f\\L^(Ba)\\‘ip\\L<^(B3ii)- 

For |x| > 2R we have that 


|xr+ 2 ®|/H=^(x)| ^ ||/||L«)(B«)[^] 5 ^(R,.)|a;r+^^ r \x + y\ 

JBr 


— n—2s 


dy 


and we remark that |j/| ^ \x\/2 (otherwise y ^ supp/). Then we use the bound 
\x + y\^ \x\ — \y\ ^ |x |/2 and we have that 


ixr+^i/5^(x)i<ii/ii«,[^°^(R„)ixr+2^ f 

J B 


I —n—2s 


dy — 


We can iterate the same method to prove that (1 + |x|"+^'*)|Il“/*^(x)| is bounded 
since D°‘f*ip{x) = f*D°‘ip{x) and D°^ip e 5s(M"). For n ^ 2s we have that 


|V;(x)| < \f{x)\\p{x)\ ^ ||/||l”(r)Ci|x 
\ip{x)\ ^ \f{x)\\p{x)\ ^ II/IIl^cr) 


2s 


for |x| ^ 1 , 
for |x| > 1 . 


Moreover, for |x| > 1 


\lp'ix)\ ^ \f{x)\\p'{x)\ + \—f{x)p{x) 


< ll/llL“(R)r4 


J/(e)(*c)< 




|<p(x)| 


Al 

<ll/llL-(R)^ + ll//(aiUHK)A<7^ 
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and for |a;| ^ 1, since / e C“(K) 

W{x)\ ^ \f{x)\\tf\x)\ + |<p(a;)| 

^||/||z,o.(K)cl|x|2-i+ci|a:p1|e/(ail l‘»(b)|2;| ^ — Clxp® 

Hence ip satisfies (2.7). Taking into account (2.23) and applying Proposition 2.4 
we have that 

f * ^{x)ip{x) dx = ^{x)ij’ix) dx = ( 27 r|a;|)“^®- 0 (x) dx, 

Jb" Jb" Jb" 


and from (2.22) we conclude that 

I f * ^{x)(f{x) dx = I {2Tr\x\)~'^^ f{x)(f{x) dx. 


□ 


The function <i> gives the representation formula for equation (—A)®m = / both 
in the distributional sense and pointwise. 

Theorem 2.8. Let f e and let u be defined as 

u{x) := ^ * f (x). 

Then u e Lj(IR") and in the distributional sense 

i-AYu = f. 

Moreover, u e C'^®+®(K") and pointwise in K” 

(-A)®u(a:) = f{x). 

Proof. From Lemma 2.6, we have that u e L^(]R"). We prove at first the statement 
for /e C®(M"). 

We notice that for (p e 5(K"), the function J'“^((27r|^)p®i^(^)) e 5s(K") and 
(27r|^)P®i^(^) satisfies the hypothesis of Lemma 2.7. Hence, by (2.20) 

< u{x), (-A) V >s= J f* <t>{x)P~^ (^(27r|^|)2"(^(5)^ (x) dx 

= f fiOpiO = f f{x)T{x) dx. 

Jfi" Jb" 

The last equality follows since / G L^(IR"), which is assured by the infinite differ¬ 
entiability of /. We conclude that u is the distributional solution of 

(-A)*w = /. 

We consider now / G C' 2 ®+'^(K”). We take a sequence of functions {fk)k e 
C®(M") such that \\fk — /||lqo(r»i) —> 0 and we consider Uk = ^ * fk- Then we 

' ' k—ftx) 

have that for any (p G 5(K") 

< Ufc, (-A)^(p >s= fk{x)(p{x)dx. 

JR^ 


By definition of fk 
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/c—»-oo 


0 . 


moreover, using (1.11) and (2.19) we have that 

<Uk-u, (-A)V [(-A)V]5,(R")ll^fc “ “IUm®") 

^ C„^s^/j[( —A) ll/fc ~ /||l‘»(R") 

We thus obtain that for any ip e 5(]R") 

< u,{-Ayp >s= f{x)p{x)dx. 

Jr" 

Hence in the distributional sense (—A)®u = / on K" for any / e Cg'*+^(M"). To 
obtain the pointwise solution, we notice that thanks to Theorem 9.3 in [13] we have 
that u e C'^®"'"®(K"). This, together with the fact that u e Lg(K”) implies that the 
fractional Laplacian of u is well defined, according to Remark 1.2. Moreover, from 
the continuity of the mapping K" B x ^ (—A)®u(a;), according to Proposition 2.1.7 
from [12], we have that (—A)'’M(a:)(/?(a;) dx is well defined. For any p e C“(M") 
we have that 

u{x){—Ayp{x) dx = f{x)p{x)dx. 

JR" JR" 

Thanks to Fubini-Tonelli’s Theorem and changing variables we obtain that for any 
(^G C”(K”) 

f{x)p{x)dx= u{x){—Ayp{x) dx = {—Ayu{x)p{x)dx. 

JR" JR" JR" 

Since both / and (—A)®u are continuous, we conclude that pointwise in R” 

= f{x). □ 


As a corollary, we have a representation formula for a (^“(R") function. 
Corollary 2.9. For any / G C®(R") there exists a function p G C®(R”) such that 

f{x) = p* $(a;), 

and p{x) = 0(1x1“"'“^®) as |x| ^ oo. 

Proof. For / G C“(R”), we define p as 

p{x) := {-Ayf{x). 

The bound established in (1.9) assures the asymptotic behavior of p, while it is not 
hard to see that p G C”(R"'). Then by using Theorem 2.8 we have that pointwise 
in R" 


p * $(x) = (—A)®/ H= $(a;) = f{x). 


□ 


2.3. The Poisson kernel. We claim that plays the role of the fractional Poisson 
kernel. Indeed, the function arises in the construction of the solution to Dirichlet 
problem with vanishing Laplacian inside the ball and a known forcing term outside 
the ball, as stated in the next Theorem 2.10. 

Theorem 2.10. Let r > 0, g e Pg(R") n (^(R") and let 

, , Pr{y,x)g{y)dy ifxeBr, 

Ug{x) := < Jr"\_b,, 

[yix) 


(2.24) 


ifxeR^\Br. 
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Then Ug is the unique pointwise continuous solution of the problem (1.1) 

J(—A)®m = 0 in Br, 

1 M = 5 in M.^\Br. 


Proof of Theorem 2.10. We see at first that Ug e Lg(K”). Take R> 2r and x e Br, 
then by using (A.22), the inequality \x — y\ > \y\ — r and for \y\ > R the bound 


(2.25) 

we have that 


\y\ 


n+2s 


< 2 


n+s 


(|i/|2-r2)®|x-i/|" 

(a:)| < f Pr{y,x)\g{y)\dy + \ Pr{y,x)\g{y)\dy 

jR>\v\>r J\v\>R 


jR>\y\>r 

^c(n,s) ^up \g{y)\ + 2'^+^c{n,s){r‘^-\x\'^y { 

yeBrt\Br. -if 

^ c(n,s) sup \g{y) \ + 2”+®c(n, s)r^® f 


\9iy)\ 


\y\>R |j/| 

\9iy)\ 


n+2s 


dy 


yeBii\B^ 


l\y\>R 


\y\u+2s 


dy. 


Since g e Lg(R”), the last integral is bounded, and so Ug is bounded in Br. It 
follows that Ug e Lg(M”), as stated. Moreover, the local C® regularity of Ug in Br 
follows from the regularity of the Poisson kernel. 

Let us fix a: e Br and prove that Ug has the s-mean value property in cc. If this 
holds, indeed. Theorem 2.2 implies that (—A)®w(a;) = 0, and given the arbitrary 
choice of x, the same is true in the whole Br. 

We claim that for any p such that 0<p<r—|a;|we have 


(2.26) Ap * Ug{x) = Ug{x). 

Let at first g be in (^“(M"). By Corollary 2.9, there exists a function ip e C”(M") 
such that 

g{y)= ^{z - y)p{z) dz 

and at infinity p{z) = (!l(|z|“"“^®). For r > 0 fixed, we write g as 

(2.27) 9{y)= f ^{t - y)(p{f) dt + f (^{z - y)<p{z) dz. 

JBr 

Using identity (A.33) we have that 

$(z - y)(p{z) dz = 

JBr 

Therefore, in (2.27) it follows that 

(2.28) 9{y) = \ <^iy - yyit) dt, 


Pr{t, z)^{y — t)dt\ (p{z) dz 


Id 

r r Pr{t,z)(p{z)dz] dt. 

V JBr J 
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where tp{t) = ip{t) + Pr{t, z)ip{z) dz. In particular, using (2.24) and (2.28) we 
have that 

= f Priy,x)( r dtj dy 

J\y\>r \ J\t\>r J 

= f V’(^) ( f Pr{y-, x)‘^{y - t) dyj dt= { i/i(t)4>(a; - t) dt 

thanks to (A.33). Furthermore, we compute 

Ap*Ug{x)= Ap{y)( - y - t) dt] dy 

J|y|>p / 

= f V'Wf f Ap{y)<i>{x - y - t) dy'] dt. 

J|i|>r VJly|>p / 

Having chosen p ^ r — |ai| we have that |a; — t| ^ |t| — \x\'^ p and from (A.31) we 
obtain 

Ap*Ug{x)= 'ijj{t)^{x — t) dt. 

J\t\>r 

Consequently Ap * Ug{x) = Ug{x), thus for g e (^“(K") the claim (2.26) is proved. 


We now prove the claim (2.26) for any forcing term g e Lg(K”) n (^(K"). In 
particular, let pk e (^“(K”) be such that gk{x) e [0,1], = 1 in Bk and pfc = 0 in 
Bk+i. Then gk := PkQ e (^“(M") and we have that gk —> g pointwise in M”, in 

k—fOD 

norm Lg(M") and uniformly on compact sets. So, for any k ^ 0 the function Ug^(x) 
has the s-mean value property in x. Precisely, for any p > 0 small independent of 
fc, 

(2.29) {Ap*Ug^){x) = Ug^{x). 

We claim that 


(2.30) 

and that for any p > 0 small 

(2.31) 


lim Ug^(x) = Ug{x) 




lim [Ap H= Ug^){x) = Ap H= ug{x). 

A ;—>00 

Let € be any arbitrarily small quantity. For k large and R > 2r, we take advantage 
of (2.25) and obtain that for x e B^ 


,{x)-ug{x)\^ \gk{y) - g{y)\Pr{y,x)dy 

JR^\Br 

^ 2”+"c(n, s)(r2 - \xfy f 




sup \gkiy) - giy)\ \ Priy,x)dy 


dy 


^ C ' 


yeBii\B, 

{n,s,r) 

JR'' 


fBR\B., 

\gk{y)-g{y)\ 




\y\ 


n+2s 


dy+ sup \gk{y) - g{y)\ ^ £ 

yeBii\Br 


by the convergence in Lg(M") norm, the uniform convergence on compact sets of gk 
to g and integrability in M.'^\Bj. of the Poisson kernel (by identity (A.22)). Hence, 
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claim (2.30) is proved. In order to prove claim (2.31), we notice that for any p > 0 
small we have that 


- Ap*Ug{x)\ ^ Ap{y)\ug^{x - y) - Ug{x - y)\dy 
4|y|>p 

^J|y|>P ^p^y)\3k{x-y) - g{x-y)\dy 

\x-y\'^r 

+ J ^3kiz) - g{z)\^ Ap{y)Pr{z,x-y)dydz 


(2.32) 


)\z\>r 

= II + l2- 


\x-y\<r 


Let R > 2p. Thanks to the bound (2.25) for \y\ > R, the convergence in Lg(M") 
norm, the uniform convergence on compact sets of g^ to g and the integrability in 
M.'^\Bp of the s-mean kernel (by identity (A. 19)) we have that for k large 

\gk{x~y) -g{x-y)\ 


h = c 


{n,sy^ f 

■^1 

(n,s,r) 


\y\>p 

\x-y\lir 


s: 2”+®c 


(^\y\2 - p2y\y\n 

\gk{x - y) - g{x - y)\ 


dy 


\y\>R 


\y\ 


n+2s 


dy 


+ ^up \gkix - y) - g{x - y)\ \ Ap{y)dy^^. 

JR>\v\>p z 


R>\v\>p 


Once more, for R> 2r and \z\ > R we use the bound (2.25) and we have that 


I 2 = 


\gk{z)-g{z)\ Ap{y)Pr{z,x-y)dydz 

J|z|>R J, 

\x—y\<r 

+ \gk{z)-g{z)\ Ap{y)Pr{z,x-y)dydz 

jR>\z\>r J, 'y'.P 


\x-y\<r 


< c(n,s) Ap{y){r‘^-\x-y\y^ 


\gk{z)-g{z)\ 


\x-^\<r 


\z\>R {\z\^-r‘^Y\z-x + y\^ 


dz dy 


+ sup 

ZE:B 


| 5 fr(z)-5(3)1 J 

(n. 5 ,r)f 

J|. 


\y\>p 

\x-y\<r 


(y) f 

Ji 


R>\z\>r 


Pr{z, X — y)dz dy 


|..ln+ 2 ^ \ 9 k{z)-g{z)\ 

\z\>R FI zeBR\Br 


since by identity (A.19) and (A.22) 


I u ^p{y) Pr{z,y - x)dzdy ^ 1. 

J. F JR>\z\>r 

\x-y\<r 


JR>\z\>r 


Therefore again by the convergence in Lg(M^) norm, the uniform convergence on 

£■ 

compact sets of gk we have that I 2 ^ In (2.32) it follows that 

lim iAp * UgA{x) = Ap* Ug{x), 
fc^GO 

thus the desired result (2.31). 
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By (2.29), (2.30) and (2.31) we have that 

Ap * Ug{x) = lim Ap * Ug^{x) = lim Ug^{x) = Ug{x), 
k —^00 k —^00 

thus Ug has the s-mean value property at x. This concludes the proof of the claim 
(2.26). 

We now prove the continuity of Ug. Of course, Ug is continuous in and in 
'MA\Br. We need to check the continuity at the boundary of Br- 

Let yo e dBr and e > 0 arbitrarily small to be fixed, (5e > 0 be such that, if 
y e Bs^{yo) then \giy) — giyo) \ < £■ We fix y arbitrarily small such that 0 < /i < ^, 
R > 2r, and x e B^ ri Bp{yo). Notice that 

— |a;p = (r + |a;|)(r — |x|) < 2r\yQ — xl < 2ry. 


From (A.22) we have that 

(2-33) |ug(a;) - g(yo)| < f \g{y) - g{yo)\Pr{y, x) dy. 

JR^\Br 

For r < \y\ < R and \y — yo\ ^ we have that \x — y\ ^ 6^ — y > Meanwhile, 
for \y\ > R we use the bound (2.25). We have that 


r 

J \v\>r 

Iv-yol^Sc 


\g{y) - giyo)\PT{y,x)dy 


/ 2 ^ r 

^c(n,s,i?)/i*(^—J 


|5(2/)I + l5(yo)| 


i?,>|y|>r 

\v-Vo\>Se 


r2^s 


dy + 2 


n+s 


I 


l'?(y)l + \9iyo)\ 




\y\ 


n+2s 


dy 


^ c{n, s, R)y" \^—c{r, R, s, g) + 2"'+®||5||ii(R„) + c{g,s,R) 
= C{n,s,R,r,g,de)y’'. 


From this and the fact that 


|y|>r- \9iy) - 9iyo)\Pr{y,x)dy \ Pr{y,x)dy = £ 

Ay-yo\<s^ 


by (A.22) and the continuity of g, we can pass to the limit in (2.33). Sending first 
y —>■ 0 and afterwards e —> 0 we obtain that 


lim {ug{x) - g{yo)) = 0, 

x^yo 


thus the continuity of Ug. 

The uniqueness of the solution follows from the Maximum Principle. Indeed, 
if one takes ui and U 2 two different continuous solutions of the Dirichlet problem, 
then It = ui — U 2 is a continuous solution to the problem 


{-Ayu{x) = 0 , 

u(x) = 0, 


in Br 
in m^\Br. 


By Theorem 3.3.3 in [3], the solution u is constant, hence null since it is continuous 
in M" and vanishing outside of Br- This concludes the proof of the Theorem. □ 
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3. The Green function for the ball 

The purpose of this section is to prove Theorems 3.1 and 3.2. Theorem 3.1 
introduces a formula for the Green function on the ball, that is more suitable for 
applications. In Theorem 3.2 the solution to the Dirichlet problem with vanishing 
data outside the ball and a given term inside the ball is built in terms of the Green 
function. We also compute the normalization constants needed in the formula of 
the Green function on the ball. 


3.1. Formula of the Green function for the ball. This subsection focuses on 
the proof of Theorem 3.1, that we recall here. 


Theorem 3.1. 

if n 2s 


Let r > 0 be fixed and let G be the funetion defined in (1.16). 

^ra(x,z) ^s-1 


G{x, z) = K(n, s)|z — x\ 


f 

Jo 


(t + 1)^ 


dt. 


Then 


where 

(3.1) 


ro(a;,z) = 




and K{n, s) is a constant depending only on n and s. 
For n = 2s, the following holds 


(3.2) 


G{x,z) = /c(^l, log 


/ r'^ — xz + ^ (r^ — x'^){r'^ 
\ r\z — ai| 



We consider the three cases n > 2s, n < 2s and n = 2s separately. 


Proof of Theorem 3.1 for n > 2s. Let x, z e Br be fixed. 

We insert the explicit formula (1.13) into definition (1.16) and obtain that 

(3.3) G{x, z) = a{n, s) (|z — — A(x, z )), 


where 


A{x,z) := r 
J\- 

(1.14) I 

A{x, z) = c{n, s) 

Jh 


Pr{y,x) 


dy. 


t\y\>r \y 

Inserting also definition (1.14) we have that 

(r2- |a;n« 


’\v\>r 


\y 


\n—2s 




r2y\y. 


■dy. 


We use the point inversion transformation that is detailed in the Appendix. 
Let X* e ]SA\Br and y* e B^ be the inversion of x, respectively y with center at z, 
defined by the relation (A.17). With this transformation, using formulas (A.18a), 
(A. 18b) and (A. 18c) we obtain that 


A{x, z) 


c{n, s)\z — X 



(|x*|2-r2)^ 

— |y*p)s|a;* — 


|n 


dy*. 


We continue the proof for n > 3. However, the results hold for n ^ 3 and 
can be proved with similar computations. We use hyperspherical coordinates with 
p > 0 and 0, 01,..., d „_3 e [O,7r],0„_2 e [0,27r] (see (A.23) in the Appendix and 
observations therein). Without loss of generality and up to rotations, we assume 
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that X* = \x*\en, so we have the identity lx* — ?/*P = + \x*\'^ — 2|x*|pcos0 (see 

Figure 2 in the Appendix for clarity). With this change of coordinates, we obtain 


^ 3 n-jz 

A(x,2) = c(n,s)|z-x|2'’-"(|x*p-r2)*27rn sin'=6>d6» 

fc=l‘^0 

r ( r sin"-^g \ 

Jo + |a:*p — 2|x*|pcos0)"-/2 j ^ 

Let T := (notice that r > 1). We have that 


/■ 


sin”-2 e 


■ de = 


1 r 

P” Jo 


sin”-2 0 


Jo (p^ + |a:*p — 2|x*|pcos0)"/2 P” Jo (t^ + 1 — 2r cos0)"/^ 
Thanks to identity (A.34) we obtain that 


d9. 


L 


sin"-^ 0 


■d0 = 


1 


0 (p^ + |x*P — 2|x*|pCOS0)"/2 pn 

1 


/■ 


sin" ^ a da 

sin"“^ada. 

Jo 


C*|"“2(|x*P — p2) Jo 

Then, using identity (A.35) and inserting the explicit value of c(n, s) given by (1-22), 
we arrive at 


(3.4) 


where 


TT 

sin(7rs), 


,* |2 


0‘ 


f 

JO 


2p 


n—1 


|a;*|"“^ Jo (r^ — p2)®(|x*p — p2) 

„(|x*|2-r2) 


dp 


In—2 


^J(x*), 


j(x*) = r 

Jo 


2p" 


0 (r^ — p2)'’(|x*p — p2) 


dp. 


Now we define the constant 


(3.5) 


k{n,s) := - 


^n—2s—l 


{l-ry-^dT 


(we compute its explicit value in Subsection 3.3). Then we have that 


J{x*) = 2k{n, s) r 

Jo 


2p" 


^(1 — T^)® ^ dr ) dp. 


0 (r2-p2)s(|a;H=|2_p2) 

We perform the change of variables t = rp and apply the Fubini-Tonelli’s theorem 
to obtain that 


J(x*) = 2k{n, s) r ^ 

Jo 

= 2fc(n, s) r 

Jo 


.n—2s —1/^2 j.2\s 


0 (r2-p2)s(|a;*|2_^2) 

2p(p2 - t2)-l 


(p^-n 


dt] 


dp 


.n—2s — l 


t (r2 — p'^)^(\x* |2 — 


dp ) dt. 
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We change variables p‘^ — = t and — t — = p to obtain 


J(a;*) = 2fc(n,s) J Q 

= 2k{n, s) r 

Jo 


^s-l 


{r^ — T — — T — t^) 


dr dt 


±rL—2s—l 


r -t (J1 _ +2 _ 


(r^ — ir — pY 


0 p^{\x* p + p — r^) 

2k{n,s) r dt, 

Jo 


dp ) dt 


where 


nt) = \ 

Jo 


(^2 _4.2 _ „\s-l 


{r^ — V' — pY 


■ dp. 


Jo p®(|a:*p + p — r^) 

Using Proposition A.11 for a = — t^ and Y = \x*Y — we have that 


m = 




sin(7rs) (|a:*p —r^)® 


Hence in J{x*), with the changes of variables = r and then — 1 = f we have 
that 


J(x*) = 2k(n, s )^-.—r(|a;' \~ — r 
sinfTTs) 


!)-. r 

Jo 


*|2_^2',-s _fTt-2s-l/i =1= |2 _ i2',s-l , 


-ty-^ dt. 


= 2k{n, s) — 


TT X 


5i«|n—2 rco / 2 -|\s —1 


= k{n, s) 


{t - 1)' 

sinjTTs) (|a;* p — r^)® jMil T-n-i 

^ \n—2 pGO 


j: 

j: 


dr 

dt. 


sin(7rs) (lx*P — r^)® JmMiM (t + 1)"/^ 
Using formula (A.18a) and definition (3.1) we have the equalities 

| 2 ;*| 2_^2 ^ j ,2 _ |^| 2 ^^ j ,2 _ |_ j ,| 2 ^ ^ ^ 


Therefore inserting J{x*) into (3.4) it follows that 

r® t 

A{x,z) = k{n,s)\z — xY^ " 

Jr 

By inserting this into (3.3) we obtain that 


ro(x,z) {t + 1)"^^ 


dt. 


G{x,z) = a{n,s)\z — xY^ — k{n,s) f 

V Jroix, 


z) (t + l)"/2 


dt I. 


Now we change the variable t = 1/r^ — 1 in definition (3.5) and obtain that 

'-CX) +S-1 

7- rjrdt= 1. 

Jo (^ + 1)^ 


(3.6) 

It follows that 


r 

k(n,s) 

Jo 


G{x^ z) = a(n, s)k{n^ s)\z — x 


2s—n 


I 


roix,z) ^s-1 


(t+1)^ 


dt. 


We set 
(3.7) 


K{n,s) := a(ri, s)k(n, s) 
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and conclude that 

G{x, z) = K{n, s)|z — 

Hence the desired result in the case n > 2s. 


i 


ro(x,z) fs-l 


(t+1)^ 


dt. 


□ 


Proof of Theorem 3.1 for n < 2s. We consider without loss of generality r = 1 (by 
rescaling, the statement of the theorem is verified in the more general case). By 
(1.13) and definition (1.16) we have that 

(3.8) G{x,z) = a(l,s)(|z - - A{x,z)), 

where 

Pi{y.x) 


A{x,z) := r 

Jr 

Using definition (1.14) we have that 
A{x,z) = c(l,s) 

Jb\(-i,i) 


(-1T) \z-y\^ - 

(l-x^)'* 


dy. 


\y 


|l-2sO,2 


(y2 _ l)s|j;-a: 


dy. 


We proceed exactly as in the case n > 2s performing the point inversion transfor¬ 
mation and we arrive at 


A{x, z) = c(l, s)|z — a; 




(x* -1)" 


(1 — y*'^Y\x* — y* 


dy* 


where |a:*| > 1. By symmetry we have that 

(3.9) A{x,z) = c{l,s)\z — x\^‘^~^{x* — l)®|x*|J(x*), 

with 


J{x*) = f 
Jo 


(1 — y* )^{x* — y* ) 


We change the variable y* = t and obtain that 


dy* 




dt. 


By the integral representation (A. 14) of the hypergeometric function, it follows that 




1 ^2r(i)r(i-s) 


1 3 


1 




ni-s) 

We use the linear transformation (A.lSd) (notice that (l/x*)^ < 1) and obtain 

that 

(3.10) 


1 3 


1 




r(|-s)r(-s) 

r(i_s)r(i-s)^ 


P{ 1; Ai S -I- 1, 


X*^ - 1 


i^-« r(f-s)r(s ) 

r(i)r(i) 


Fl --S, 1-s, 1-s, 


x*"^ - 1 


The first hypergeometric function obtained in the sum (3.10) is transformed ac¬ 
cording to (A. 15c) as 




x*"^ - 1 


= -,s,s-l-l,l — a;* 
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For s + 1 > s > 0 and 11 — x* | < 1 the convergence conditions are fulfilled for the 
integral representation (A. 14) of the hypergeometric function. Therefore we may 
write 

,r(s + i) 


1; Aj S + Ij 


- 1 


= X- 


L 


r(s) Jo 


1/2 


dt. 


On the other hand, for the second hypergeometric function obtained in identity 

(3.10), we use transformations (A. 15b) and (A. 15c) and arrive at 


F 


1 


s, 1 - s, 1 - s. 


x*^-l 




1 


s, 0,1 — s, 1 — a;* 


ai=^|A 0,-,1-s, 


- 1 


We use the Gauss expansion (A.11) with a = 0, 5=^, c=l — s and w = 

(we notice that 0>c—a — 6>—lfors> 1/2 and |i(;| < 1, thus the series is 
convergent). Since a = 0, all the terms of the series vanish, except for fc = 0. Hence 
we obtain that 




1 


= 1 


and therefore 


Consequently 


I rr-*^ 

Fl - — s, 1 — s, 1 — s, 


X*" - 1 


= x*~^"lx*l 


1 / r(i)r(-s)r(s + i) r(s)r(i-s) 

1^1 V r(i-s)r(s) Jo (i + (3;.2_i)/)V2 

We recall that c(l, s) = (r(s)r(l — s)) ^ and we define the constant 

(3.11) fc(l, s) := c(l,s)- 


^r(i)r(-s)r(s + i) 


r(i-s)r(s) 


We insert J(x*) into (3.9) and have that 


A(a;, z) = k(l, s)\z — x\ 


-'I 


(a 


lyF 


0 (1 + (cc*^ - l)t) 


1/2 


dt + \z — x\' 


With the change of variables (a:*^ — l)t = t we obtain that 


A(x, z) = fc(l, s)|z — x\ 


-I 




{t + 1) 


Y dt + I z — XI 


Inserting this into (3.8) and noticing that x* — 1 = ro(x, z) it follows that 


G(x, z) = —a(l, s)k(l, s)|z — x\ 


croix,z) ^-s-l 

-1 I ^ 

(t + 1) 


I 


Y dt. 


We call 

(3.12) /c(l, s) = —a(l, s)fc(l, s) 

and conclude the proof of Theorem 3.1 for n < 2s. 


□ 
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Proof of Theorem 3.1 for n = 2s. Without loss of generality, we assume r = 1. We 
insert the explicit formulas (1-13) and (1-14) into definition (1.16). Moreover, we 

use the explicit values of the constant a^l, from (1.21) and c^l, from (1.22). 


We obtain that 
(3.13) G{x,z) =--\og\x - z\ +\ 


Let 


r \og\y-z\J ^ ^ 
J\v\>i V 2^ 1 

A{x,z):=\ \og\y-z\A- 

J\y\^i V 2 


dy 


’- 2 / 1 ' 


y'^ - 1 \x - y\' 

We perform the change of variables v = ^i^ice 1 — ^ 0, we have that 

|u| ^ 1. We set w := and observe that |w| ^ 1. It follows that 


A{x,z)=\ flog^^ + log|z-x| 
V 


u — x\ 


dv 

7r^ 


We use identity (A.32) and since |?r;| 5= 1 and |a;| ^ 1 we obtain that 
A{x, z) = Trlog ^|w| + {w^ — 1)^^^^ + Trlog \x — z\ 

= 7rlog(l — zx + '\/(l — a;^)(l — z^)). 
Inserting this into (3.13) we obtain that 


' I — zx + -^(1 — x'^){l — z^) 
\x-z\ 

This completes the proof of Theorem 3.1 for n = 2s. 


G{x,z) = ^log 

TT 


□ 


3.2. Representation formula for the fractional Poisson equation. This sub¬ 
section is dedicated to the proof of Theorem 3.2, which we recall here. 

Theorem 3.2. Let r > 0, he C'^®+'^(ilr) n G{Br) and let 

, . Ky)G{x,y)dy ifxeBr, 

u{x) := < Jb, 

[o ifxeW^\Br. 

Then u is the unique pointwise continuous solution of the problem (1.2) 

{ (—A)^m = h in Br, 

= 0 in M.^\Br. 

Proof of Theorem 3.2. We identify h with its extension, namely we con¬ 

sider h e (72®+®(M’^) with B^ Cl supp/i such that h = h on Br. Then, by definition 
(1.16) we have that in Br 

u(x) = h(z)G(x,z)dz 

JBr 

= h(z)^(z — x)dz — h{z)( ^{y — z)Pr{y,x)dy\dz 

JBr JBr \JTSL^\Br J 

= /iH=$(x)- Pr{y,x)[h*^){y)dy. 
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Let 

g{x) := h* $(a:) for any x e K". 

As we have seen in Theorem 2.8, g e n C(M"). Let for any x e ^ 

u{x) = Vo{x) - vi{x), 

where Vo{x) = g{x) in K" and 


vi{x) = 


I 

JR' 




Pr{y,x)g{y)dy lixsBr, 


^g{x) if X e 

Then for x e Br, thanks to Theorems 2.8 and 2.10 

(—A)®m(x) = h{x) — 0 = /i(x), 

hence u is solution (1.2). Also, Theorem 2.8 and Theorem 2.10 assure the continuity 
of u in K". 

The uniqueness of the solution follows from the simple application of the Maxi¬ 
mum Principle for the fractional Laplacian (see Theorem 3.3.3 in [3]). □ 

3.3. Computation of the normalization constants. This subsection is dedi¬ 
cated to the computation of the constant K,{n, s) in Theorem 3.1. 

We first compute the constant k{n,s) given in (3.5) for n > 2s. We claim that 


(3.14) 


k{n, s) = 


r(f) 


r(t-s)r(s)- 


Indeed, using definition (3.5) and taking the change of variable = t we have 
that 


1 


k{n, s) 


= 2 


J 

Jo 


^n—2s — l 


(l-r^) 


.2\s-l 




dt = \ 12 
0 


s —1 


(1 —t)® dt. 


We use identities (A.9) and (A. 10) to obtain that 




t 2 


-S-1 


(1 — ty ^ dt = 


r(f - s)r(s) 
r(f) ^ 


which is exactly the result. 

We now prove Theorem 3.3, namely we compute the constant K{n, s) encountered 
in the formula of the Green function G. 

Theorem 3.3. The constant K,{n,s) introduced in identity (1.17) is 

r(f) 


c(n, s) = 


22'*7rtr2(s) 




for n 7 ^ 2s, 
for n = 2s. 


Proof of Theorem 3.3. For n > 2s, we insert the values of a{n,s) from (1.20) and 
of k{n,s) from (3.14) into (3.7) and we obtain that 


K(n,s) = a{n, s)k{n, s) = 


r(t) 


22«7r2r2(s) 
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For n < 2s, we recall definitions (3.11), (3.12) and (1.20), we use identities (A.5), 
(A.7) and (A.3) relative to the Gamma function and obtain that 


t(l,s) = —a(l, s)A:(l, s) = 


(-s)r(-s) 


1 


1 


22«r(s) r(i-s)r(s) 22«r2(s)' 


On the other hand, we recall that ^^(^1, as we have seen in the proof of 

Theorem 3.1 for n = 2s. This concludes the proof of Theorem 3.3. □ 


We make now a remark on the constants C{n,s) defined by (1.6) and c(n,s) 
defined in (1.22). These two constants are both used in various works in the defi¬ 
nition of the fractional Laplacian, but as we have seen in the course of this paper, 
they arise for different normalization purposes. The constant C(n, s) as defined by 
[4] is consistent with the Fourier expression of the fractional Laplacian, meanwhile 
c(n, s) as introduced in [10] is used to normalize the Poisson kernel (and the s-mean 
kernel), and is consistent with the constants used for the fundamental solution and 
the Green function. We first introduce the direct computation of the fractional 
Laplacian of a particular function. Namely: 

Lemma 3.4. Let u(x) = (1 — |a;p) + . Then in Bi 

{-Aru{x) = C{n,s)‘^B{s,l-s), 

where B is the Beta function defined in (A.8). 

The more general case (more precisely, for the function u{x) = (1 — |a:p)(( for 
any p > —1) was proved in [6, 5]. With small modifications with respect to the 
general case, the proof of Lemma 3.4 can be also found in Section 3.6 in [3]. 


Theorem 3.5. The constant C{n,s) introduced in (1.6) is given by 

2^"sr (§ -t s) 


(3.15) 


C{n, s) = 


TT 2 F(1 — s) 

Proof of Theorem 3.5. By Lemma 3.4 we have that in Bi 




{-AYu{x) = C{n,s)—B{l - s,s). 
We use Theorem 3.2 and for n ^ 2s, we obtain that 


u{x)= r C{n,s)‘^B{l-s,s)G{x,y)dy 


- 


= C{n,s)^B{l-s 


\x - y\ 

We compute the latter identity in zero and have that 


r / rro{x,y) fs-1 \ 


1 — l-y I • ^ 

(3.16) 1 = C(n,s)^B(l - s,s)K(n,s) dt^ dy. 
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We compute the double integral, by using Fubini-Tonelli’s theorem 


JBi \Jo 


|yP t 


s-1 




■\y-it 


t 


s-l 


(t + 1) 


■w dt\ dp 




UJn 

Ys 


0 (^ + 1 ) 

00 ^s-l 


r 


UJ. 


(t+1)' 




By inserting this, the value of «;(n, s) from Theorem 3.3 and the measure of the 
(n — 1)-dimensional unit sphere a;„ = (27r"’/^)/r(n/2) into (3.16) and using (A.10) 
we obtain that 

22^sr(f + s) 

C{n,s) = —IT— - 

TT 2 i (1 — 

For n = 2s we have that (—A)^u(x) = C (1,1/2) tt. Thanks to Theorem 3.2 


u{x) = C ( 1, i ) TT 


J G{x,y)dy. 


Using formula (3.2) and computing u at zero, we obtain that 


1 = C 1,:, 


a)i: 


log 


1 + y/l - 

\y\ 


dy = ttC ( 1, - ) . 


Hence C (1,1/2) = I/tt and this concludes the proof of the Theorem. 


□ 


Appendix A. Appendix 

A.l. The Gamma, Beta and hypergeometric functions. We recall here a few 
notions on the special functions Gamma, Beta and hypergeometric (see [1], Chap¬ 
ters 6 and 15 for details). 


Gamma function. The Gamma function is defined for x > 0 as (see [1], 
Chapter 6): 

pOO 

(A.l) r(x) := dt. 

Jo 

The Gamma function has an unique continuation to the whole K except at the 
negative integers, by means of Euler’s infinite product. We have that r(l) = 


r( 2 ) = 

1 and r(l/2) = y'fr. We also recall the next useful identities: 

(A.2) 

r(n -1- 1) = n! 

for any n e N, 

(A.3) 

r(x + 1) = xT{x) 

for any x > 0, 

(A.4) 

r(l /2 -t x) v^ 2 i- 2 “ 
r( 2 x) r(x) 

for any x > 0, 

(A.5) 

r(s)r(i-s) 

sin(7rs) 

for s e (0,1), 

(A.6) 

r(i /2 - s)r(i/ 2 -p s) = 

cos(7rs) 

for s e (0,1), 

(A.7) 

T{1-s) = {-s)T{-s) 

for s e (0,1). 
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Beta function. The Beta function can be represented as an integral (see [1], 
Section 6.2), namely for x,y>Q 


rOO j.X — 1 

(A.8) (1+ !)»» ■“ 

and equivalently 

(A.9) B{x,y) = r — t)y~^ dt. 

Jo 

Furthermore, in relation to the Gamma function we have the identity 


(A.IO) 


B{x,y) 


r(cc)r(i/) _ 

r(a; + y) ' 


Hypergeometric functions. There are several representations for the hyper¬ 
geometric function (see [1], Chapter 15). We recall the ones useful for our own 
purposes. 


(1) Gauss series 


(A.ll) 


00 

F{a,b,c,w) = 2 

A;=0 


{a)k{b)k w’^ 
{c)k k\ ’ 


where {q)k is the Pochhammer symbol defined by: 


(A.12) 




1 for fc = 0, 

q{q -f 1) • • • (g -I- /c — 1) for fc > 0. 


The interval of convergence of the series is |w| ^ 1. The Gauss series, on its 
interval of convergence, diverges when c — a—b ^ —1, is absolutely convergent when 
c — a — b > 0 and is conditionally convergent when |t(;| < 1 and —l<c—a—6^0. 
Also, the series is not defined when c is a negative integer —m, provided a or 5 is 
a positive integer n and n < m. 

Some useful elementary computations are 


(A.13a) 

(A.13b) 


F{a, b,b,w) = {1 — w) 


1 


1 


A( a, - -I- a, -, w ) = 


(1 -I- w) -f (1 — w)~ 


2 ’ 2 ’ 

(2) Integral representation 

(A.14) ^^(0,6,0,?!;) := - t)“"'’"i(l - wf)"^ df. 

The integral is convergent (thus F is defined) when c> b > 0 and |u>| < 1. 


(3) Linear transformation formulas 
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From the integral representation (A. 14), the following transformations can be 
deduced. 

(A.15a) 


F{a,b,c,w) = {1 — wY “ ^F{c—a,c—b,c,w), 


(A.15b) 

(A.15c) 


{l-w)-^F 

(a,c b,c, ^ 

V w — 1/ 

1 

1 

1—1 

b,c a,c, 1 

^ w — 1/ 


(A.lSd) 


r(c)r(c — a — b) 
r(c — a)r(c — b) 

+ (1 


F(a, b, a + b — c + 1,1 — w) 

a.c — b.c — a ■ 




r(a)r(i) 


6+1,1 — w), 


when 0 < w < 1. 


A.2. Point inversion transformations. The purpose of this appendix is to recall 
some basic geometric features of the point inversion, related to the so-called Kelvin 
transformation. 

Let r > 0 to be fixed. 


Definition A.l. Let xq b Br he a, fixed point. The inversion with center xq is a 
point transformation that maps an arbitrary point y e IR”\{xo} to the point Ji.xgljj) 
such that the points y, xq, lie on one line, xq separates y and K 2 ;(,(y) and 

(A.16) K^Yv) ■= ^0 - - ^^^{y-xo). 

\y-xor 

This is a bijective map from IR”\{a;o} onto itself. Of course, Ka;j,(Ka;Q (a;)) = X. 
When this does not generate any confusion, we will use the notation y* := Uxoiu) 
and X* := Ka;Q(a:) to denote the inversion of y and x respectively, with center at 
a;o. 


Remark A.2. It is not hard to see, from definition (A.16), that 
(A.17) \y* - xo\\y - xo\ = Y - Ixol"^■ 


Proposition A.3. Let Xq e Br he a fixed point, and x* and y* be the inversion 
of X e K"'\{xo} respectively y e K"\{xo} with center at Xq. Then: 


a) points on the sphere dBr are mapped into points on the same sphere, 

b) points outside the sphere dB^ are mapped into points inside the sphere, 
(A.18a) 


f'l \y-xo? ^ 1 

^ (r2 - |xo|2)(r2 - |2/|2) |y=f|2_r2’ 

(A.18b) 

d) 

\y - xoY \y* - xoY' 

(A.18c) 

e) \y*-x*\ = {r^-\xo\^)r-^^^ 

\y - xo\\x - 


a:o| ’ 
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The Kelvin point inversion transformation is well known (see, for instance, the 
Appendix in [10]) and elementary geometrical considerations can be used to prove 
this lemma. We give here a sketch of the proof. 

Sketch of the proof. A simple way to prove claims a) to c) is to consider the first 
triangle in Figure 1. 



Figure 1. Inversion of x,y with center at xq 

We denote b := \OY\ = \y\, b* := |OF*| = |i/*|, a := [AoKj = \y — xq\ and 
P := |AoK*| = \y* — Xoj. Let OH be the perpendicular from O onto YY* . We 
apply the Pythagorean Theorem in the three triangles AOYH, AOHXq, AOHY* , 
add the equation (A.17) and by solving the system, one gets that 

^*2 _ 

a 

From this, claims a) to c) follow after elementary computations. 

In order to prove d), without loss of generality, one can consider the point inver¬ 
sion of radius one with center at zero y* = —y/\y\'^ and take its derivative. Since 
the point inversion transformation is invariant under rotation, we can assume that 
y = jyjei and the desired result plainly follows. 

To prove e), see the Appendix in [10], or consider the second triangle in Figure 
1. We denote a := [AqA] = \x — xq], b := \XqY\ = \y — xq], a := \XY\ = \x — y\ 
and p := jA*y*j = jx* — y*\. Let YHi and Y*H2 be perpendiculars from Y, 
respectively Y* onto the segment AA*. By applying the Pythagorean Theorem 
in the four triangles AXqYHi, AXYHi, AXqY*H2, AX*Y*H2^ adding relation 
(A.17) and using that YHi is parallel to Y*H2, one gets after solving the system 
that 

_ (r^ - lxol^)a 

^ ab 

which is the desired result. □ 

A.3. Some useful integral identities. We present here a few detailed computa¬ 
tions related to the functions $, and Pr and some other useful integral identities. 
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Lemma A.4. For any r > 0 
(A.19) 


K{y) dy = 1 . 


Proof. Using (1.12) and passing to polar coordinates we have that 


Ar{y)dy = c(n, s) 

Jm.' 

fO 

= c(n, s)a;„ 

Jr 


r.2s 


B, ilyl"^ -r^y\y\' 


dy 


p{p^ — r^)' 


dp, 


where lT{nft) is the measure of the (n — l)-dimensional unit sphere. 

We change the variable z = (p/u)^ — 1 and have that 


(A.20) 


J 

JR-' 


'\Br 


Ar{y) dy = 


c(ji, s) 


-w, 


I 


0 iz+l)z‘ 


■ dz. 


We apply identities (A. 8 ) and (A.10) and use identity (A.5) to obtain that 

1 

JO {z+Oz'^ 


(A.21) 


r 

JO 


dz = r(i — s)r(s) = 


TT 

sinTTs 


Using the definition (1.22) of c(n, s) it follows that -^r{y) dy = 1, as desired. 

□ 


Lemma A.5. For any r > 0 and any x e Br 

(A.22) r Pr{y,x)dy = l. 

Proof. From the definition (1.14) of Pr we have that 



We make the proof for n > 3. However, the results hold for n ^ 3. We 
change variables using the hyperspherical coordinates with radius p > 0 and angles 
6 », 6 »i,..., 0 n -3 e [ 0 , tt] , 6 I „_2 e [0,27r] 

yi =p sin 0 sin 0i... sin Ons sin 0„_2 
2/2 =P sin 0 sin 0i... sin 0„_3 cos 0n-2 
(A.23) 2/3 =psin 6 »sindi... cos 6»„_3 


yn =pcos0. 

The Jacobian of the transformation is given by sin"“^ 0sin”“^ 0i... sin 0 „_ 3 . 
We only remark that for n = 3 the usual spherical coordinates can be used 2/1 = 
p sin 0 sin 01 , j /2 = psin 0 cos 0 i and 2/3 = pcosd, while for n = 2 and n = 1 similar 
computations can be performed. 

Without loss of generality and up to rotations, we assume that x = |x|e„ to 
obtain the identity \x — y\'^ = p^ + \x\^ — 2|a;|pcos0 (see Figure 2 for clarity). With 
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this change of coordinates, we obtain 
Pr{y,x)dy 

n—3 


=c(n, s)(r^ 


2)*27r 


It — O pTT pOU r 

sin*^6»d6» 

-^0 Jr Jo 


p"-isin”-2 6»d6»dp 


0 (/o^ — + |a;p — 2/9 |x| cos0)"/2 


We do the substitution r = r/\x\ and p = p/\x\ but still use r and p for simplicity 


6,1 



Figure 2. 


and we remark that now p > 1 and r > 1. We obtain that 
(A.24) 

Pr{y,x)dy 


n —3 ^TT ^00 f r 

(n,s)(T^-ir2Tn{sin'=0d0[ (7^(1 


sin""2 0d£/ 


Let 


*(P) := f 

Jo 


(p2 _ j.2y (^p2 I _ 2pcos6*)^/^ 


dp. 


0 — 2pcos0 + 1)”/^ 


We claim that, for p > 1 

(A.25) i{p) = ^2 _ 1) edO. 

To prove this, we use the following change of coordinates 

sind sin Of 

(A.26) 

We have that 
(A.27) 


yj p^ — 2p cos 9 + 1 P 


de= 1 - 


cos a 


V? 


da. 


sm a 
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To see this, one takes the derivative of the relation (A.26) 
(A.28) 


(pcos0 — l)(p — COS0) cos a , 
dO = - da 


— 2pcos6 + 1)2 p 

and obtains with some manipulations of (A.26) that 

cos a\/p^ — sin^ a p(pcos6> — l)(p — cos0) 


p^ — sin^ a — cos a 
Now by changing variables we obtain that 


— 2/9 cos 6* + 1) 2 


*(P)= f 

Jo 


sin"-2 6» 


0 (p^ — 2p cos 9 + l)"/2 


de 


„n-2 


f 

JO 


sin" a da 


(V - si 


sm^ a — cos a)-\/p^ — sin^ a 


1 sin" ^ q;(V p"^ — sin^ a + cos a) da 


„n-2 


L 


1 


— 1)V p^ — sin^ 


P" - 1) V Jo 


f 


sin" ^ a da + 


L 


sin" ^ a cos a 


■\Jp^ — sin^ a 


da 


By symmetry 


therefore 


I 


^ sin" ^ a cos a 
\/p^ — sin^ a 


da = 0, 




1 


f 

Jo 


pn-2(p2 _ 

We substitute this into (A.24) and obtain that 

n—2 


sin" a da. 




-\B. 

We claim that 
(A.29) 


n—2 


76 — Z p- 

nj 

k=l ->0 


sin^ 6d9 = 


n/2 


r(n/2)- 


To prove this, we integrate by parts and obtain that 


pTT nTi nTi 

9 d9 = [k — 1) sin^“^ P dP — (fc — 1) si 

Jo Jo Jo 


4 = sin'^PdP ={k- 

Jo 

which implies that 


sin'^PdP, 


Thus we have 


Ik = 


Ik = " 


k-l 


r sin'^ ^ PdP = 

Jo k 


fc-lP-3 1 

k-lk-3 2, 
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with Jo = TT and Ji = 2, and the claim (A.29) follows after elementary computations. 
And so 

7r"/2 


J 

JR' 


"\Br. 


Pr{y,x)dy = c{n,s){r - 1)'* 


J 


2p 


r(n/2) (p 2 _ r 2 )«(p 2 _ 1 ) 


dp. 


We change variable = -sr and obtain 


r 7r"/2 r® 

= c(n,s) 

Jh"\b, r(n/2) Jo 


0 + 1) 


dz. 


We use (A.21) and the value of c{n, s) from (1.22) and obtain that 

Pr{y,x)dy=l. 

This completes the proof of Lemma A.22. 

Lemma A. 6. For any r > 0 and any x B Br 
(A.30) -- f 


□ 


:(b,s) 

JL 


B, - |yp)1a; - yl" 


= 1 . 


Proof. Let y* be the inversion of y with center at x (notice that |y*| > r). Then 
by using (A. 18a) and (A. 18b) we obtain that 


L 


dy 


J 

Jr' 


r2 — |a:|2 


dy* 


\B, I |2/*P -r2 / \x-y*\^’ 


Jb, - y|"-2^ 

From identity (A.22) the desired result immediately follows. 
Lemma A.7. For any r > 0 and any x G W^\Br 


□ 


(A.31) 


Ar{y)^{x-y)dy = ^{x). 

Jr’*\b^ 


Proof. We prove the claim for n 2s. We insert definitions (1.12) and (1.13) and 
obtain that 


Ar{y)<^{x-y)dy = r^‘'c{n,s)a{n,s)\ 
jR"\Br JR' 


1 


AS. 


— r2)®|y|"|x — y\^ 


■dy. 


Let X* and y* be the inversion of x, respectively y with center at 0. Using identities 
(A.17), (A.18a) (A.18c) and (A.18b) we obtain that 


Ar{y)<^{x - y) dy = 

Jr"\b. 

From (A.30) it follows that 


c(n, s)a(n, s) 


In—2s 


I 


dy* 


Br \x* — y*|"“ 2 s ^j .2 _ 


Ar{y)^{x - y) dy = 

Jr"\b. 


a(n, s) 


In—2s ’ 


and thus the desired result. 

We now prove the claim for n = 2s, assuming r = 1. We have that 

r A.(y)$(x-y)dy = ^ r } ^J^^ ^Ldy. 
Jr"\b. tt J|y|>i vy - i|j/I 
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We perform the change of variables v = ^, with jul ^ 1. We set w := hence 
I re I ^ 1. Then we have that 


f Ar{y)‘^>{x - y) dy = ^ \ ("log + log |a;|") 


dv 


VT^ 


We use the following result (see [2], page 549) 

dv 


(A.32) r 

Jl 


lr.aU,_nl_:2_ = i if|a|^l 

I„|s 5 i ® Vl - [Trlogdal + (a^ - 1)^/^) - 7rlog2, if |a| > 1. 


We thus obtain 

Ar{y)<^{x - y) dy =--log\x\, 

which concludes the proof of the Lemma. □ 

Lemma A.8. For any r > 0, let xq e Sr he a hxed point. For any x e W^\Br 


(A.33) 


Pr(y,xo)^(x - y) dy = ^(x - xo). 


Proof. We prove the claim for n A 2s. We have that 


Pr{y,xo)^{x-y)dy = c{n,s)a{n,s) \ 

JIR 


'\Br 


(r^ - |aiop)dai - dy 

(|y|2-r2)dy-a;o|"' 


Let X* and y* be the inversion of x, respectively y with center at xq. From (A. 17), 
(A. 18a) (A. 18b) and (A. 18c) we have that 


Pr{y,xo)^{x - y)dy 


-\B. 

= c{n, s)a{n, s)~ 


|x* -xo|”"2s 


£ 


dy* 


(r2 — |xop)"“^* Jsr ~ \y*\‘^)^\y* — X* |"“2s 

Using (A.30), we obtain that 

r Pr{y,xo)<^{x-y)dy = 
which concludes the proof for n ^ 2s. 


|n-2s ’ 


\X - Xo| 

We now prove the claim for n = 2s, assuming r = 1. We have that 


r Pr{y,xo)<^{x-y)dy = ^ I J 

J|y|>l V 


1 - Xq^ log 1^ - x\ 
j/2 - 1 \y-xo\ 


dy. 


We perform the change of variables v = ■, noticing that jul ^ 1. We set 


w := hence Iwl ^ 1. Then we have that 

X — Xq ’ ' ' 


\ 

JR'- 


Pr{y,Xo)^{x -y)dy = 


"\Br 

We use identity (A.32) and obtain 




1 \v-w\ , . 

log]-1 +log|x-Xo| 

k) - xo 


dv 


Vl — u2 


which concludes the proof. 


r Pr(y,a;o)4’(x-j/)dy =--log|x-xo| 


□ 
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We emphasize here two computations that we used in the proof of Lemma A.5, 
namely identities (A.25) and (A.29). 


Proposition A. 9. For any t > 1 

r sin"-" e 

Jo (r^ — 2 tcos0 + 1)"/^ 

Proposition A. 10. 

n-2 

(A.35) TT 


de = 


nn—2 {^2 


n—2 pTY 

nf 

7._1 JO 


sin^ 6 d9 = 


(r" - 1) 


rt/2 


f 

Jo 


sin" a da. 


k = l 


T{n/2y 


In the next Proposition we introduce yet another integral identity. 
Proposition A. 11. Let a,/3 e K such that < 1- Then 

{a — xY~^ TT (a +/3)®-^ 


r 

Jo 


dx = 


x®(/3 + x) sin(7rs) 

Proof. We change the variable x = at and obtain that 

r = A f' ~ 

Jo xYP + x) /?jo V Pj 

We use the integral definition (A. 14) of the hypergeometric function for a = 1, 
&=1 — s, c=l and w = —j (since |t| < 1 the integral is convergent) and we 
obtain that 


J' 

Jo 




-1 


dz = 


r(s)r(i-s) 

r(i) 


F( l,l-s,l,-- j. 


Now we use the linear transformation (A. 15c) and compute 

Ffl-s,0,l, 


a 




a + fi 

~T~ 


a 


a + /3 


We use the Gauss expansion in (A.11) and notice that for k > 0, all the terms of 
the sum vanish. We are left with only with the term k = 0 and obtain that 


F 1 - s,0,l, 


a 


a + P 


= 1 . 


Furthermore, F(1 — s)r(s) = (as in identity (A. 5)) and it follows that 


f 


(a — xY ^ , TT (a + jdY ^ 
dx = —-T-, 


YP + x) 


□ 


sinTTs 


We explicitly compute here another integral that was used in our computations, 
namely : 

Proposition A. 12. For any s e (0,1/2] we have that 


(A.36) 


r<X> 

sin t dt = —cos{7rs)r{2s — 1). 
Jo 
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Proof. We have that 


(A.37) 


f 


sint dt =-Im I dt. 


\ 


±2s—2^—it 


We consider the closed path rip = d^([0, p\ x [0, p\) n -Bp(O)^ . We take the contour 
integral dz, and let 7 p = dBp(Q) n ([0,p] x [0,p]) (the boundary of 

the quarter of the circle). By Cauchy’s Theorem, the contour integral is 0 (there 
are no poles inside flp), therefore 


Jo hp 


z'^^-^e-^dz- 


JO 


dt = 0. 


Integrating along 7 p, by using polar coordinates z = pe*® and then the change of 
variables cos 9 = t we have that 


r 


C7r/2 


C^/2 

z^^-^e-^dz 

= 



e-P^°^^d9 

J^p 


Jo 


Jo 


= P 


I 


1 g-pt 




dt 




+ cp^®“^ 

p1/2 

e-^‘ dt 


Jl/2 


Jo 


= ^ -^2s-2(g-p/2 _ 


Hence 


z2«-2e-" dz = 0 


lim 

and we are left only with the integrals along the real and the imaginary axis, namely 

rCO pOO 

t2s-2g-i ^ ■2s-l t2s-2g-it 

Jo Jo 


Here the left hand side returns the Gamma function according to definition (A.l). 
We compute = (cos(7r/2)+i sin(7r/2)) = sin(7rs)+jcos(7rs) and in (A.37) 

we obtain that 

J sintdt = —r(2s — l)Im^sin(7rs) + icos(7rs)^ = — cos(7rs)r(2s — 1). 

This concludes the proof of the Proposition. □ 
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